Quantum critical point shifts under superconductivity: 
the pnictides and the cuprates 
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FIG. 1: The SDW and nematic critical points shifts from 
Fermi surface effects. The tuning parameters r, s are for the 
SDW and the nematic phase transitions. The critical points 
in the metal are at r^,s^, and under superconductivity, these 
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We compare the position of an ordering transition in a metal to that in a superconductor. For 
the spin density wave (SDW) transition, we find that the quantum critical point shifts by order |A|, 
where A is pairing amplitude, so that the region of SDW order is smaller in the superconductor 
than in the metal. This shift is larger than the ~ |Ap shift predicted by theories of competing 
orders which ignore Fermi surface effects. For Ising-nematic order, the shift from Fermi surface 
effects remains of order |A| . We discuss implications of these results for the phase diagrams of 
the cuprates and the pnictides. We conclude that recent observations imply that the Ising-nematic 
order is tied to the square of the SDW order in the pnictides, but not in the cuprates. 



I. INTRODUCTION 

The interplay between spin density wave (SDW) or- 
dering and superconductivity clearly plays a central role 
in the physics of a variety of quasi two-dimensional cor- 
related electron materials. This is evident from recent 
studies of the phase diagram of the ferro-pnictidesir— and 
the '115' family of heavy- fermion compounds^^. In the 
cuprates, it has been argued that d-wave superconduc- 
tivity is induced by SDW fluctuations in a metal^, and 
this has been the starting point for numerous studies of 
the complex phase diagram-'^^'^^. In all these materials, 
there is a regime of co-existence between SDW ordering 
and superconductivity, and this opens the way to a study 
of the 'competition' between these ordersi^: this compe- 
tition can be tuned by an applied magnetic field, as has 
been studied in a number of revealing experiments^""— 
on the LSCO and YBCO series of superconductors. 

This paper will discuss a question that arises naturally 
in the study of such competing orders2ii^. We consider, 
first, the 'parent' quantum critical point as that associ- 
ated with the onset of SDW order, (^, in a metal. To ac- 
cess this point we have to suppress superconductivity in 
some manner, say by the application of a magnetic field. 
This parent critical point will occur at a value r^ of some 
tuning parameter r, which could be the carrier concentra- 
tion or the applied pressure. We define r so that r < r" is 
the SDW phase with (i^) 7^ 0; see Fig. [TJ The value of r° 
is clearly material specific, and will depend upon numer- 
ous microscopic details. Then, we turn our attention to 
the onset of SDW order within the superconductor (SC); 
we characterize the latter by a gap amplitude A, and 
denote the critical value of r by r^ . The essence of the 
picture of competing orders is that the onset of super- 
conductivity should shrink the region of SDW order, and 
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towards the ordered phases. 



the 115 compounds^. In the pnictides, a "backbending" 
of the onset of SDW order upon entering the SC phase, 
consistent with the idea of r" — 7-^ > 0. 

Let us begin by computing the shift r° — r^ in mean- 
field Landau theory. The simplest free energy of the SDW 
and SC order parameters has the form^: 

/:-(r-r,")<^2^7/(^2^' + f|A|2+u|A|4 + K(^2|A|2. (1) 

Here k > is the phenomenological parameter which 
controls the competition between the order parameters. 
Examining the onset of a phase with {ip) 7^ in the 
superconductor with A 7^ 0, we conclude immediately 
from Eq. (P that 
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hence 
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We will be interested here in particular 



in the magnitude of the shift r^ — r^- We will see that 
the shift is dominated by low energy physics, and so has 



a universal character. This shift 



played a central 



role in the phase diagrams presented in Refs. i22il23l , and 
applied to the cuprates. Recent work has shown that 
similar phase diagrams also apply to the pnictides^i^ and 



Such a shift was a key feature of the theory in Ref. [231. 

The primary focus of the previous work was in the 
lower field region, where the superconductivity is well- 
formed, and A is large. Here it is appropriate to treat 
the superconductivity in a mean-field manner, and ig- 
nore pairing fiuctuations, while treating spin fluctuations 
more carefully. The present paper turns the focus to 
higher fields, where eventually superconductivity is lost. 
Here, clearly. Landau theory cannot be expected to apply 
to the superconducting order. Moreover, we expect the 
Fermi surface of the electrons to be revealed, and a more 
careful treatment of the electronic degrees of freedom is 



called for. One of the primary results of our paper will be 
that the Landau theory result in Eq. ^ breaks down for 
small A, and in particular in the limit |A| —!■ 0. This is 
a consequence of the crucial importance of Fermi surface 
physics in determining the position of the SDW transi- 
tion at T = 0. Instead, we will show from the physics 
of the "hot spots" on the Fermi surface that the shift is 
larger, with 



CIAI 



(3) 



For the competing order picture to hold, we require that 
C > 0. Somewhat surprisingly, we will find that our 
results for C are not transparently positive definite. Dif- 
ferent regions of the Fermi surface contribute opposing 
signs, so that determining the final sign of C becomes a 
delicate computation. In particular C will depend upon 
the vicinity of 'hot spots' on the Fermi surface, which are 
special points connected by the SDW ordering wavevec- 
tor. We will find that the immediate vicinity of the hot 
spots contributes a positive sign to C, while farther re- 
gions contribute a negative sign. Thus the primary com- 
petition between SDW and superconductivity happens 
at the hot spots, while other regions of the Fermi surface 
which survive the onset of SDW order continue to yield 
an attraction between SDW and superconductivity. For 
the case where the two hotspots connected by the SDW 
ordering wavevector are equivalent under a lattice sym- 
metry operation (i.e. they have the same pairing gap 
and the same magnitude of the Fermi velocity) , we will 
find that distinct contributions to C exactly compensate 
each other, so that C = 0. However, in the case that the 
two spots are not crystallographically equivalent (which 
is the generic situation in both the cuprates and the pnic- 
tides), we will show that C > 0. A positive C is indicated 
in Fig. m 

We had considered the shift in SDW ordering due 
to superconductivity in a previous workSi. However, 
in that work, the metallic and superconducting states 
were not Fermi liquids and BCS states respectively, but 
rather fractionalized states known as 'algebraic charge 
liquids—""—. In this case, we found that the competition 
between SDW and superconductivity was robust, and al- 
ways yielded a shrinking in the size of the SDW region. 
We will not consider such exotic states here, but work 
entirely within the framework of Fermi liquid theory, in 
which the onset of superconductivity leads to a tradi- 
tional BCS superconductor. In this context the interplay 
between SDW and SC in a Fermi liquid is conveniently 
encapsulated in the 'spin-fermion' model^. We will find 
the same qualitative shift in the SDW critical point as 
found earlier^^, and the estimate in Eq. ([3]). 

We will also generalize our methods to analyze the 
shifts in the quantum critical points of other orderings be- 
tween the metallic and superconducting phases. Specifi- 
cally, we will consider charge density wave (CDW) order 
and Ising-nematic order, 77. We will find that the CDW 
shift initially appears to be formally similar to Eq. ([3]), 



but the coefficient C is found to be exactly zero; terms 
higher order in A do indicate competition with super- 
conductivity, but the CDW critical point shift is much 
smaller than the SDW's For Ising-nematic order, 77, we 
will also find that the Fermi surface result is similar to 
that in Landau theory, as in Eq. ([5]). This smaller shift 
is also illustrated in Fig. [TJ The weaker effect upon 77 
is due to its reduced sensitivity to the gap opened by 
superconductivity on the Fermi surface. 

We will present a detailed discussion of the implica- 
tion of these results for the pnictide and cuprate phase 
diagrams to Section IVII However, let us highlight here 
an important inference that will follow from our com- 
putations. We will argue that the experimental phase 
diagrams imply that Ising-nematic ordering and SDW 
ordering are independent instabilities of the Fermi sur- 
face for the cuprates. In contrast, for the pnictides, our 
conclusion will be that the SDW ordering is the primary 
Fermi surface instability, and the Ising-nematic ordering 
is a secondary response to the square of the SDW order 
parameter. 

This paper is structured as follows. We will begin in 
Section [H] by introducing our starting point, the spin- 
fermion model with the SDW fluctuation. Within the 
spin-fermion theory, we show that the SDW fluctuation 

induces the d wave pairing for the cuprate and s^ for 

the pnictides instead usual s wave pairings in Section iHll 
Assuming the d or s^ wave pairings in each case, we ex- 
tend the spin-fermion theory into the theory with pairing 
and other possible orders in Section ITVl In section |Vl we 
show the quantum critical point shifts toward the or- 
dered phase, which explicitly shows the competition be- 
tween superconductivity and the SDW phase. Section IVll 
presents our conclusions. 



II. THE SPIN-FERMION MODEL 

We will study the system with SDW quantum phase 
transition in two dimensional system. The main ingre- 
dients of the spin-fermion model are Fermi surfaces and 
the SDW order parameter. Let us first consider generic 
microscopic Hamiltonians for the cuprates and the pnic- 
tides. 

Hcu = y^ec(fc)4ca + Hsdw 
k 

Hpe = Yl '^c{k)cica + Y^ €d{k)dlda + H^dw (4) 
fc k 

Here, we consider "minimal Fermi surfaces" for both ma- 
terials, where c correspond to the one band electron for 
the cuprates around the F point, and for the pnictides, 
c, d correspond to the hole and electron bands centered at 
(0, 0) and (±7r/ao, 0). All terms containing the SDW op- 
erators are in Hgdw In Fig- HI we illustrate typical hole- 
doping cuprates large Fermi surface and pnictides' two 
band Fermi surfaces. To see general features, we consider 



incommensurate ordering wave- vectors for cupratei 



,20 



Qi 



27r 



1 



Q2 



27r 



1 1 
2' 2 



(5) 



For the pnictides, we represent the ordering wavevector 
as Q = (tt/oojO) explicitly in the figure, but there is an- 
other hot spot with the ordering vector, Q = (0,7r/ao). 
From now on, we set the lattice constant as a fundamen- 
tal unit as usual. As it is well-known in the literature,— 
the incommensurate SDW fluctuation is decribed by two 
complex vector wave functions. 
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In the figure, two distinct hot spots are represented by 
filled and empty circles in the cuprate Fermi surface. The 
filled one is farther from the nodal point than the empty 
one. Note that the incommensurate SDW fluctuation 
links a filled circle with a empty circle. If we consider one 
special case, the commensurate SDW fluctuation, two 
kinds of hot spots become identical, and we only need 
one real 0(3) field to describe the SDW fluctuation as 
usual. 



Because the hot spots mainly contribute to the tran- 
sition in the low energy theory, we write the electron 
annihilation operator as combination of the hot spot con- 
tinuum fields as follows. For the cuprate. 



, j = (l,---,4,I,...,4). 



(7) 



The /, g fields represent the nearer and farther fields from 
the nodal points, so the SDW fluctuation mixes the / and 
g fields. The commensurate limit means / and g fields 
become identical. For the pnictides case, two bands can 
be described by hot spot fields as follows. 
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As in the cuprates, the SDW fluctuation mixes / and 
g particles in this notation. We note that the electron 
band does not have the shape of circle, which makes the 
SDW possible, and there is no symmetry that guaran- 
tees the identity of the two bands. For convenience, we 
focus on the commensurate case in the remaining of this 
section, and the next two sections for describing spin- 
fermion models. But, later in the Sec. |Vl we will come 
back to the general incommensurate cases and consider 
the critical point shifts in general. 

The spin-fermion model with commensurate SDW sim- 
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FIG. 2: Hot spots near SDW transition. In the upper and 
lower panel, large Fermi surface for the cuprates and two band 
structure for the pnictides are shown. The ordering vectors 
for the cuprates are Qi = |^(i - ^, i), Q2 = |f (i § - ^)), 

and the pnictides have Q = (7r,0). In the upper panel, 
the filled and empty circles have different distances from the 
node, which means the gap magnitudes are different. Two 
ordering wave vectors are represented by the dashed and thin 
lines in the upper panel and arrowed lines are in the bot- 
tom. For the pnictides, the electron band is distorted because 
there is no symmetry that guarantees identicalness of the hole 
and the electron band, which can also induces different gap 
functions.— 1^ In both panels, every hot spot is numbered and 
the bar notation is used for the negative. 



ply becomes 
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where M,- 



nected by Q. The first and second lines describe the dy- 
namics of the spin and fermion sectors. The third line is 
the "Yukawa" coupling term. Note that we explicitly in- 
clude the second derivative kinetic term for the fermions 
in the theory, which becomes irrelevant if we only focus 
on the SDW phase transition. Such term is not neces- 
sary for describing the SDW phase transition with non- 




the pnictides, the pairing and its vertex correction are 



FIG. 3: Pairing vertex. The line with arrow is fermion and 
the wavy line is for the SDW fluctuation. 



coUinear Fermi velocities only, but the existence plays an 
important role in extending the theory to the one with 
pairing and nematic orders. Note that the final form of 
the spin-fermion model is exactly the same in both the 
cuprates and the pnictides even though the microscopic 
band structure and the ordering vectors are completely 
different. This means the physics for the SDW transition 
is universal and we can focus on one case and apply the 
result to the other case. 

As usual, the effective action for the SDW phase tran- 
sition is given by integrating out the fermions and expand 
with order parameters. 
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This Hertz- type theory is well-known and it describes the 
SDW fluctuation with dynamical critical exponent z = 2 
at least in the zeroth order. In this paper, we only focus 
on the critical point shifts rather than critical properties 
of the transition itself.— 
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Note that the relative sign between hot spots does not 
change, which is the main characteristic of the s wave 
pairing. As we can see, the vertex becomes irrelevant 
in the low energy limit in the RG sense. Therefore, the 
SDW fluctuation cannot mediate usual s wave pairing for 
the cuprates and s++ for the pnictides. 

However, in the d wave channel of the cuprates^ and 

the s-i channel^iiH of the pnictides, the relative sign 

between the hotspots changes, given the pairing symme- 
tries. Such relative sign changes allow the pairing chan- 
nel's enhancement. The pairing and its vertex correction 
are 
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Clearly, the alternative sign change induces enhancement 

of the d, s.| wave pairings in the low energy limit. 

Therefore, the d and s^ pairing is natural under the 

SDW fluctuations rather than usual s and s_|_+ pairings 
in conventional theory. In the next sections, we assume 
the existence of the pairings in each case and incorporate 
them into the spin-fermion model in a manner consistent 
with symmetry. Moreover, by symmetry consideration, 
we also introduce other possible order parameters such 
as a nematic order and charge density wave and extend 
our theory to incorporate them. 

Notice that the vertex correction in the Eq. ([T2l) is log- 
arithmically divergent if we have finite correlation length. 
Such behavior is a well-known signature of the conven- 
tional BCS theory. However, in the quantum critical re- 
gion, the pairing boson is softened and the quantum fluc- 
tuation becomes important. The nature of this quantum 
critical pairing has been discussed in Refs. [llll28l 



III. PAIRING INSTABILITIES 



Within the spin-fermion model, we can address pairing 
problems naturally. If we consider the SDW fluctuation 
as a pairing boson, then we need to investigate plausi- 
bility of the pairing instabilities by the SDW. The basic 
idea is following. If we assume there is infinitesimal pair- 
ing, then the pairing becomes enhanced or suppressed by 
the integrating out higher energy-momentum contribu- 
tions depending on the possibility of the pairing channel. 
In the Fig. [3] such a pairing vertex is illustrated. Note 
that the fermions with opposite momentums are paired, 
so the participating fermions in the pairing is not the 
same as ones in the SDW in general. If we consider the 
s wave channel in the cuprates and the s_|_+ channel for 



IV. MICROSCOPIC SYMMETRY AND 
EFFECTIVE THEORY 

To extend the spin-fermion model to the one with pair- 
ing terms and nematic order parameter, let us consider 
microscopic symmetries and thier transformations be- 
cause the square lattice symmetry should be respected 
in the low energy theory. Hereafter, we analyze the sym- 
metry in terms of the cuprate problem unless otherwise 
stated. It is easy to extend it to the pnictides case. Due 
to the dx2-y2 wave property, the pairing term's rotation 
and reflection needs additional factors. In Appendix [^ 
we show the explicit transformation properties of fields 
and bilinear terms to avoid notation ambiguity. Nambu 
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TABLE I: Symmetry transformations of the spinor fields un- 
der square lattice symmetry operations. Tx,y: translation by 
one lattice spacing along the x, y direction; R-^/2'- 90° rotation 
about a lattice site (a; — >■ j/, y — >■ —a;); Ixy'- reflection about 
the X — y axis {x ^ y,y ^- x); T: time-reversal, defined as a 
symmetry (similar to parity) of the imaginary time path inte- 
gral and the conjugate fields are transformed to ^^ —5- ^Vy. 
Note that such a T operation is not anti-linear. Also, the 
notation, ^i^b = (^J j,)^ is used for convenience. 



spinors for particles are defined in a usual way: 
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(13) 



In Table m we summarize the transformation rules of the 
spinor fields. 

Among various combinations of bilinear terms, the fol- 
lowing operators are of interest. 

-^L^'^3,a + *la^'^4,a (W) 

+ *,L^"*3,a + *i,a^"*4,a (15) 

-4.^'*3,a + <,r^*4,a (16) 

Op ^ +*L^^*3,a + 4,r^*2,a (17) 

The operator, Oa, is invariant under all the lattice sym- 
metries. This is just d wave pairing term's low energy 
expression. Therefore, inserting the operator, Oa, in the 
original Lagrangian is certainly allowed by the symme- 
try consideration. For the future convenience, let us write 
down the fcrmion Hamiltonian with pairing explicitly. 
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ei{k) = Vi ■ k + a v'jk'^. Note that A is a given constant 
here. 

The 0^,r; operators are also interesting, and they are 
transformed as follows: 
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Therefore, O^^ri operators have nematic ordering symme- 
tries. The difference between two operators are that Oc^ 
is from the pairing channel and O,, is from the density 
channel. 

The final operator Op describes a charge density wave 
order parameter, which has horizontal ordering wavevec- 
tor in this case. Such ordering was considered a candidate 
of the "pseudogap" phase and the high energy 4ao order- 
ing pattern in the cuprates^. Notice that the hot spots 
are in general not linked by (7r/4, 0), but linked by the 
ordering vector, Qcdw — ^/2 ~ fc/3, which is consistent 
with the observations of the Hudson and collaborators'^^. 
For a thorough consideration of the charge density wave, 
we need to investigate "hot-spots" for the charge density 
wave order and start from the beginning. But because 
the calculations are identical in both cases, we consider 
the charge ordering within the present SDW theory. 

The original theory can be extended by introducing 
the other order parameters from the above microscopic 
consideration. The total theory is 
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(20) 



where Ho is the spatial representation of the Eq. ^TE\\ 
and mi^4 = l,?7i2.3 = —1 ■ Note that the C terms are 
higher order in terms of the pairings. 

Before going further, let us remark the meaning of the 
above consideration. In this extended theory, we have 
found the two different nematic order channels. One or- 
der parameter, 77, couples to the density of fermions. The 
other one, ^, couples to the pairing, which suggests they 




one ordering wavevector,say 1, is 



FIG. 4: Susceptibility Feynman diagram. The wavy line is 
for the order parameters like the SDW or the nematic order. 
The vertex matrix elements depend on the order parameters. 



have different charges. Many previous works focus on the 
pairing channel nematic order with the nodal fermions 
of the d^2_y2 superconductivityi^^— The nematic or- 
der naturally induces different pairings, A^; ^ Ay, which 
corresponds to the condensation of ( in this case. The 
original nodal fermions are gapped and the nodes be- 
come shifted depending on the sign of the gap function. 
Within the nodal fermion theory the density channel 
is not allowed by the square lattice symmetry. In our 
spin-fermion model, the density channel is surely allowed 
and can see the effect under the small superconductiv- 
ity. Note that the pairing channel is higher order in the 
pairing amplitude A, which we will assume it is small 
here. 



V. CRITICAL POINT SHIFTS 

In this section, to study influence of the superconduc- 
tivity on quantum critical points in general, we come back 
to general cases, including the incommensurate cases for 
the cuprates. We can easily generalize the previous dis- 
cussion to the incommensurate case: the changes are 
mainly in the fermions' spectra, which have different 
velocities and gap functions between two SDW linked 
points, and the existence of the complex two vector fields 
for the SDW. We will see effects of superconductivity on 
the quantum critical points by evaluating the diagrams 
in Fig. |4] for various order parameters with and without 
superconductivity. 
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Due to the symmetry it is enough to choose one of the 
linking hot-spots. For example, the SDW fluctuation be- 
tween kf^i and kg 2 is 
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The two gap functions arise from the two bands in pnic- 
tides, and from the distance differences from the nodal 
point in the cuprates. To parametrize gap difference, let 
us introduce one parameter, a as Aj- = A(l±a); thus a 
characterizes the distinction between the two hot spots, 
which become crystallographically equivalent a = 0. For 
the cases under consideration here, a = for a commen- 
surate (tt, tt) wavevector for the SDW ordering for the 
cuprate case, while a 7^ in all other cases. As we showed 
in the Eq. (|10p . the susceptibility function contributing 
to the critical point and the relative critical point shifts 
with and without superconductivity is defined as 



Xa = Xo - (Xo - Xa) 



(22) 



A positive br implies the critical point shifts to shrink the 
SDW region leading to competition between the SDW 
and the superconductivity, while the negative sign im- 
plies attraction between the SDW and the superconduc- 
tivity. 



A. Spin density wave 



To investigate the effect of the superconductivity on 
the SDW critical point, let us focus on the SDW order 
near the criticality with and without superconductivity. 
The interaction between the SDW order and fermions 
naturally affects the critical point as we saw in the Eq. 
pO| . If we allow the superconductivity, the main effect 
of the superconductivity is gapping out the fermion sur- 
faces, so the contribution to the critical point has to be 
changed. The amount of the change can be obtained by 
evaluating the susceptibility with and without the pair- 
ing. In each case, the critical point is affected with the 
loop contribution and we can evaluate them with the 
Hamiltonian, Eq. P^ . The total susceptibility for the 



As we can see, the fermion loop formula, Xa' ^^ ^^~ 
vergent in the ultraviolet limit without the curvature 
term. But the critical point shift, which was defined as 
difference between different pairing magnitudes, is well- 
defined. This is because the susceptibility function con- 
tains two independent momentum components unless the 
two Fermi velocities are parallel, which requires hot spots 
are the same as the nodes. We exclude such special case 
in this paper. Therefore, the curvature term, a, can be 
dropped, and we assume |A|a <C 1. So, the dispersion 



relation ei{k) — vl ■ k will be used for evaluation. 
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where the Fermi velocities are defined as 

vi =-uji (008(6*1), sin(0i)) , V2 = w/2(cos(6'2),sin(6'2) 

Therefore, the angle dependence and the relative gap 
function determine the critical point shift. The angle de- 
pendence on the two Fermi velocities indicates that more 
parallel velocities implies a larger critical point shift. So 
the perpendicular Fermi velocities of SDW participating 
fermions have the smallest critical point shift. As we 
can see, the sine function dependence indicates that the 
collinear Fermi velocities are dangerous in our calcula- 
tion. In this case, we need to keep a from the start, and 
the functional behavior becomes the same as the nematic 
case. The coefficient, C{a), is 



CA(a) = -T-^ / dq^dqy 






X 1 + 



Vg2 + (l-a)2^g2 + (1 + a) 



(24) 



As pointed out in Section HI the sign of CA(a) is not im- 
mediately evident from Eq. (l24l) : the small q region near 
the hot spot contributes a positive sign, while that from 
large q has a negative sign. A numerical evaluation of 
Eq. (P^ shows that the result is indeed non-negative for 
all a, and the result is shown in Fig. [5] Ci^ipt) increases 
monotonically with increasing a. 

A curious feature of Fig. [5] is that Ca(0) = 0: it van- 
ishes for the case of equivalent hot spots. However, it 
is not at all evident from the integral in Eq. ([M)) that it 
equals zero for a = 0. This is more easily proved from the 
original expression in Eq. (|21|) . by reversing the order of 
frequency and momentum integration. Evaluating first 
the momentum integration in Eq. (I21[) (after linearizing 
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FIG. 5: The SDW critical point shift with the relative gap 
parameter a, Ci\(a). 



the dispersion about the Fermi surfaces) we find 



2w/iw/2|sin(6'i - 6*2)1 

u? + A+A_ 



Ai,/c^2 + A2 



(25) 



It is now evident that the critical point shift vanishes for 
A+ = A_, for then the above result becomes indepen- 
dent of the value of the gap. 

If we consider the state slightly off the criticality. e.g. 
by considering finite temperature, then the low energy 
physics is governed by the temperature. In such a case, 
the shift becomes an analytic function of the supercon- 
ducting gap, which means quadratic gap functions scaled 
with another energy scale, instead of the linear gap. 
Therefore, we can understand the non-analytic behavior 
of the linear gap as a property of the quantum critical 
points. 



B. Charge density wave 



The CDW ordering operator, Op, with a specific order- 
ing wave vector, Qcdw = fc/2 — ^/s, was introduced in 
the previous section. The CDW has different character- 
istics from the SDW. For example, instead of the spin 
dependent vertex, it has density type operators linking 
two hot spots, and also the linked hot-spots' gap func- 
tions have the same pairing sign and magnitude. Com- 
bining all of these, the susceptibility for the CDW with 



superconductivity is 
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(26) 



The CDW critical point shift is 
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(27) 



where the Fermi velocities are defined as 



V2 = w/2(cos(6i2),sin(6'2)) , vz = '!;/2(cos(6'3), sin(6'3), 

where the angle between ^2,3 is (^2 — 6*3), and the veloci- 
ties are not collinear in general. Interestingly, the calcu- 
lation itself is quite similar to the SDW's formally. How- 
ever, a key difference is that the participating fermions 
have the exactly same gap functions and Fermi velocities 
for the CDW, and this is guaranteed by symmetry. Con- 
sequently, the prefactor in Eq. (P71) is CA(a — 0) which 
is zero. Thus the shift in the CDW critical point vanishes 
at this order. 

We can estimates higher order in A by evaluating 
Eq. (|27p with a finite momentum cutoff A. This intro- 
duces dependence of the result on A/A, which is evalu- 
ated in Appendix IbI We find a net competing effect, but 
this is formally higher order in A and so parametrically 
smaller than the SDW case. 



point correction can be evaluated from the fermion loop 
calculation as before, and it is 
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(28) 



Note that there is a crucial difference in this integration 
compared to the previous critical point shifts. Because 
the nematic order parameter consists of particle and hole 
with same Fermi velocity, the pathology of collinear dis- 
persions are always present in the nematic phase transi- 
tion. Therefore, we need to keep the curvature term, a, 
and then we have the susceptibility as 
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(29) 



where Dq is the constant density of state. The lower cut- 
off is from the dispersion relation Sk — v-k + av'^k'^. Note 
that this integration is well-defined in both ultra-violet 
and infra-red regions. Even in the SDW the colinear 
Fermi velocity hot spots whose cases are excluded in this 
paper also suffer similar problems. 



The critical point shift of the nematic ordering, then. 
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(30) 



C. Nematic order 

With the extension of the spin-fermion theory, we can 
consider the nematic order parameter within the theory. 
As we mentioned before, there are two channels for the 
nematic order, but for the case of a small pairing gap, it 
is enough to consider the density channels. The critical 



where Cjj — -■ Note that the nematic critical point shift 
does not contain the linear gap behavior like the SDW. 
Instead it starts from the second order and it is analytic 
in terms of the gap function. Such a term describes usual 
competing term of the Landau-Ginzburg theory, as was 
discussed in Section HI Parametrically the nematic order- 
ing is more stable than the spin density wave under the 
"weak" superconductivity. 
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FIG. 6: Cuprate phase diagram adapted from Ref. |21||22| . 
Here x is hole doping, Xm. is the position of the SDW crit- 
ical point in the metal, and Xa is the SDW critical point in 
the insulator. The shift between Xm and Xs represents the 
consequence of Eq. ([3]). The regions with 'fluctuating Fermi 
pockets' have renormalized classical thermal fluctuations of 
SDW order. Ising-nematic ordering is expected for T < T* , 
and consistent with Fig. [T] this regime is not sensitive to the 
onset of superconductivity. 



VI. CONCLUSIONS 

Our main results, summarized in Fig. [TJ have a natu- 
ral application to the physics of the cuprates. For SDW 
ordering, we have shown that there is a large shift in the 
quantum critical point due to the onset of superconduc- 
tivity, represented by Eq. ([3]). In contrast, for the Ising- 
nematic order, when considered as an independent order 
parameter, there is a significantly smaller shift, of order 
that expected in Landau theory in Eq. ^ . These results 
provide a natural basis for the phase diagram proposed in 
earlier work2ii2^, which we reproduce here in Fig. [S] The 
large shift in the SDW ordering between the metal and 
the superconductor is represented by the arrow from Xm 
to Xs ■ We assume that the Ising-nematic ordering has an 
onset around Xm, and this is barely shifted by the onset 
of superconductivity, as implied by Fig.[T] Consequently, 
long-range Ising-nematic ordering can survive for x > x^^ 
as is indicated in Fig.[Bl These results are consistent with 
recent observations of Ising-nematic ordering^^— in the 
hole-doped cuprates. 

In an ap plied magnetic field, as was discussed in 
Refs. I22l23l . the point Xg eventually merges with x^ , so 
that the SDW transition in the high-field normal state 
takes place at x = Xm- Given this, we expect that the 
Ising nematic transition will also merge (or become very 
close to) with the SDW transition at high fields. 

A notable feature of Fig. |6] is the "back-bending" of 
the crossover line bounding the region where there are 
'renormalized classical' fluctuations of local SDW order: 



100 
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FIG. 7: Phase diagram for Ba[Fei_2,Coa;]2As2 from Refs. [^, 
6. The back-bending of the SDW ordering transition in the 
superconducting phase is similar to that of Tsdw in Fig. [S] 
Here, rather than renormalized classical SDW fluctuations, 
we have true long-range order indicated by 'AFM'. The Ising- 
nematic order is present in the phase labeled 'Ort' and absent 
in that labeled 'Tet'. Note that, unlike the cuprates, the Ising- 
nematic transition follows the SDW ordering transition in the 
superconducting state too. 



this region is bounded by the line labeled T* in the nor- 
mal state, and by the line labeled Tsdw in the supercon- 
ducting state. As we have argued^ii^, this is a natural 
consequence of the shift of the in SDW quantum critical 
point from Xm to x^- 

Turning to the pnictides, we note that the back- 
bending of the SDW ordering has been clearly seen in 
recent experiments, as shown in Fig. [71 Here, because 
the stronger 3-dimensionality of the crystal structure and 
the commensurate wavevector, the region of renormal- 
ized classical SDW fluctuations becomes a region of true 
long-range order, and so is more easily detected by neu- 
tron scattering. However, Fig. [7] differs from the phase 
diagram in Fig. |6] in one important aspect: note that 
the Ising-nematic transition in Fig. [7] closely tracks the 
SDW transition in both the normal and superconducting 
states, rather than separating from it in the supercon- 
ducting state as in Fig. [6] This means that the shift in 
the nematic ordering transition due to superconductivity 
is not significantly smaller than that of the SDW tran- 
sition. This is in conflict with the situation outlined in 
Fig. [1] where the nematic transition hardly shifts relative 
to the SDW transition. 

This difference between our computations and the 
pnictide phase diagram in Fig. [7] implies that the Ising- 
nematic transition in the pnictides is not an independent 
instability associated with the electrons near the Fermi 
surface. For if it were, our computations show that it 
would barely notice the onset of superconductivity. We 
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now argue that the phase diagram in Fig. [7] can be un- 
derstood if we assume that the Ising-nematic ordering is 
primarily induced by its coupling to the square of the 
SDW order— 1^. Thus, in addition to the couphng of the 
nematic order, 77, to the fermions in Eq. (|20p . we need to 
add its couphng to the SDW order: 



-r;-ip 



>'V {0l 



4) 



(31) 



where tpx {0y) is the SDW ordering at wavevector (7r,0) 
((0, tt)). Then a correction of order | A| to the SDW fluc- 
tuations from the onset of superconductivity, wiU feed 
into a similar correction to the Ising-nematic fluctua- 
tions via a perturbation theory in A. Thus our conclu- 
sion is that A is the dominant coupling which induces 
Ising-nematic order in the pnictides. A similar conclu- 
sion appears to have been reached recently by Kimber 
et al0 based upon their analysis of the STM observa- 
tions of Chuang et al.^'^ . In contrast, for the cuprates, 
the influence of A appears significantly weaker. 
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TABLE II: Symmetry transformations of the hot-spot fields 
under square lattice symmetry operations. 
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Appendix A: Symmetry 

In this section, we set the notation for the symme- 
try transformation of the square lattice. Mean field 



Hamiltonian with the d, 
Hmf — 



x^ —y^ 



needs the pairing term as 



XI *y (^laCja + h.c) - t^y , Ci,Ci,a 



^II4a^ 



<y> 



J2KA'"''^^-^jb) + h.c. 



(Al) 



factor i, which makes the d wave pairing term invariant. 
After writing the lattice fields with continuum field, we 
can obtain the transformation in Table. |TT1 It is worth- 
while to mention that at low energy or long-wavelength 
scale the fields at hot spots can be treated as independent 
fields. 



Time reversal symmetry is obtained with low energy 
fields instead of the lattice fields, so we do not consider 
it here. See the caption of the Table. ID 



The first line describes usual hopping terms on the square 
lattice, which gives the Fermi surface. We exclude special 
'nesting' type Fermi surfaces and assume there are points 
linked by the spin density wave ordering vector. 
We start with lattice field transformations. 



^x,y 

R-n/2 



I. 



xy 



Ca{x) -!■ c'^ix') = Ca{x) 
Caix) -!■ c'^{x') = iCa{x) 
Caix) -> C^{x') = iCa{x) 



(A2) 



The rotation and reflection transformation attaches the 



*,; 






, 1 = 1,2,3,4 



(A3) 



i,b/ 



After introducing the Nambu spinors, bilinear spinors' 
transformations can be done easily. 



Physical quantities are described with bilinear terms 
such as density and pairing interactions. Below several 
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important bilinear terms are listed up to constants. 

*la^"*l,a = llaha-lljla 

*I,a^'*l,a = flafla+flha 

*L^"*l,a = £abiflaflb + hbfla) 

*l,a^''^a6e6c*2,c = -{fL^abht + fL^abflb) 

'^l.a'r^^abebc^2,c = i{fL^abhb - lla^abflb) 

*l,a'r°'^ab*l,fc == fla^abflb + fl^<fabflb 

^{,a'^''^ab'^l.b ^ fla^abflb - fi^^abflb 

*I,aT"^CTab*l,fc = ifLfL^bc + hdflbead)^ab 

*l,a^'''?ab*l,fc = {-i){llaflc^bc- hdhbead)Sab 



(A4) 



Appendix B: The CDW critical point shift 
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In Fig. [SI we illustrate the function, F/^{9), with different 
cutoffs. As we can see the larger cutoffs make the smaller 
Ca(0), even though it is positive. 



As we saw in the above, the critical point shift for the 
CDW is determined by the function. 
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are getting smaller with larger cutoffs. 



34, 8190 (1986). 

^° Ar. Abanov, A. V. Chubukov, and J. Schmalian, Advances 
in Physics 52, 119 (2003). 

" M. A. Metlitski and S. Sachdev, arXiv: 1005. 1288 

^2 S. Sachdev and S.-C. Zhang, Science 295, 452 (2002). 

^^ S. Katano, M. Sato, K. Yamada, T. Suzuki, and T. Fukase, 
Phys. Rev. B 62, R14677 (2000). 

^^ B. Lake, H. M. R0nnow, N. B. Christensen, G. Aep- 
pli, K. Lefmann, D. F. McMorrow, P. Vorderwisch, 
P. Smeibidl, N. Mangkorntong, T. Sasagawa, M. Nohara, 
H. Takagi, and T. E. Mason, Nature 415, 299 (2002). 

^^ J. M. Tranquada, C. H. Lee, K. Yamada, Y. S. Lee, 
L. P. Regnault, and H. M. R0nnow, Phys. Rev. B 69, 
174507 (2004). 

^^ B. Khaykovich, S. Wakimoto, R. J. Birgeneau, M. A. Kast- 
ner, Y. S. Lee, P. Smeibidl, P. Vorderwisch, and K. Ya- 
mada, Phys. Rev. B 71, 220508 (2005). 

^■^ J. Chang, Ch. Niedermayer, R. Gilardi, N. B. Chris- 
tensen, H. M. R0nnow, D. F. McMorrow, M. Ay, J. Stahn, 
O. Sobolev, A. Hiess, S. Pailhes, C. Baines, N. Momono, 
M. Oda, M. Ido, and J. Mesot, Phys. Rev. B 78, 104525 
(2008). 

^* J. Chang, N. B. Christensen, Ch. Niedermayer, K. Lef- 
mann, H. M. R0nnow, D. F. McMorrow, A. Schneidewind, 
P. Link, A. Hiess, M. Boehm, R. Mottl, S. Pailhes, 
N. Momono, M. Oda, M. Ido, and J. Mesot, Phys. Rev. 
Lett. 102, 177006 (2009). 

^^ D. Haug, V. Hinkov, A. Suchaneck, D. S. Inosov, 
N. B. Christensen, A. Ivanov, T. Keller, C. T. Lin, and 
B. Keimer, arXiv: 1008.4298). 

^° D. Haug, V. Hinkov, Y. Sidis, P. Bourges, N. B. Chris- 
tensen, D. S. Inosov, Ch. Niedermayer, P. Bourges, 
Y. Sidis, J. T. Park, A. Ivanov, C. T. Lin, J. Mesot, and 
B. Keimer, 

^1 E. G. Moon and S. Sachdev, Phys. Rev. B 80, 035117 
(2009). 

22 S. Sachdev, Physica Status Solidi B 247, 537 (2010) 
(arXiv:0907.0008| . 



2^ E. Demler, S. Sachdev and Y. Zhang, Phys. Rev. Lett. 

87, 067202 (2001); Y. Zhang, E. Demler, and S. Sachdev, 

Phys. Rev. B. 66, 094501 (2002). 
2'' R. K. Kaul, A. Kolezhuk, M. Levin, S. Sachdev, and T. 

Senthil, Phys. Rev. B 75 , 235122 (2007). 
25 R. K. Kaul, Y. B. Kim, S. Sachdev, and T. Senthil, Nature 

Physics 4, 28 (2008). 
2® V. Galitski and S. Sachdev, Phys. Rev. B 79, 134512 

(2009). 
2'^ S. Sachdev, M. A. Metlitski, Y. Qi, and C. Xu, Phys. Rev. 

B 80, 155129 (2009). 
2* Ar. Abanov, A.V. Chubukov, and J. Schmalian, Adv. 

Phys. 52, 119 (2003). 
2^ R. Khasanov, D. V. Evtushinsky, A. Amato, H. H. Klauss, 

H. Luetkens, Ch. Niedermayer, B. Biichner, G. L. Sun, 

C. T. Lin, J. T. Park, D. S. Inosov, and V. Hinkov, Phys. 

Rev. Lett 102, 187005 (2009). 
^° P. Szabo, Z. Pribulova, G. Pristas, S. L. Bud'ko, P. C. Can- 
field, and P. Samuely, Phys. Rev. B 79, 012503 (2009). 
^^ I. I. Mazin, D. J. Singh, M. D. Johannes, and M.-H. Du, 

Phys. Rev. Lett. 101, 057003 (2008). 
^2 Kangjun Seo, B. A. Bernevig, and Jiangping Hu, Phys. 

Rev. Lett. 101, 206404 (2008). 
^^ Y. Kohsaka, C. Taylor, K. Fujita, A. Schmidt, C. Lupien, 

T. Hanaguri, M. Azuma, M. Takano, H. Eisaki, H. Takagi, 

S. Uchida, and J. C. Davis, Science 315, 1380 (2007). 
^^ W. D. Wise, M. C. Boyer, Kamalesh Chatterjee, 

Takeshi Kondo, T. Takeuchi, H. Ikuta, Yayu Wang, and 

E. W. Hudson, Nature Physics 4, 696 (2008). 
^^ M. Vojta, Y. Zhang, and S. Sachdev, Phys. Rev. Lett. 85, 

4940 (2000); Erratum 100, 089904(E) (2008). 
^^ M. Vojta, Y. Zhang, and S. Sachdev, Int. J. Mod. Phys. B 

14, 3719 (2000). 
^'^ M. D. Prato, A. Pelissetto, and E. Vicari, Phys. Rev. B. 

74, 144507 (2006). 
^* A. Pelissetto, S. Sachdev, and E. Vicari, Phys. Rev. Lett. 

101, 027005 (2008). 
^3 E. A. Kim, M. J. Lawler, P. Oreto, S. Sachdev, E. Fradkin, 

and S. A. Kivelson, Phys. Rev. B 77, 184514 (2008). 
"^ Y. Huh and S. Sachdev, Phys. Rev. B 78, 064512 (2008). 
"^ V. Hinkov, D. Haug, B. Fauque, P. Bourges, Y. Sidis, A. 

Ivanov, C. Bernhard, C. T. Lin, and B. Keimer, Science 

319, 597 (2008). 
*2 R. Daou, J. Chang, D. LeBoeuf, O. Cyr-Choiniere, F. Lal- 

iberte, N. Doiron-Leyraud, B. J. Ramshaw, R. Liang, D. 

A. Bonn, W. N. Hardy, and L. Taillefer, Nature 463, 519 

(2010). 
*^ M. J. Lawler, K. Fujita, Jhinhwan Lee, A. R. Schmidt, 

Y. Kohsaka, Chung Koo Kim, H. Eisaki, S. Uchida, 

J. C. Davis, J. P. Sethna, and Eun-Ah Kim, Nature 466, 

347 (2010). 
'"' C. Fang, H. Yao, W. Tsai, J. Hu, S. A. Kivelson, Phys. 

Rev. B 77, 224509 (2008). 
"'^ C. Xu, M. MiiUer, and S. Sachdev, Phys. Rev. B. 78, 

020501 (R) (2008). 
""^ S. A. J. Kimber, D. N. Argyriou, and I. I. Mazin, 

arXiv: 1005. 1761 (2010). 
"■^ T.-M. Chuang, M. P. Allan, Jinho Lee, Yang Xie, Ni Ni, S. 

L. Bud'ko, G. S. Boebinger, P. C. Canfield and J. C. Davis, 



Science, 327, 181 (2010). 



